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We propose a method to achieve significant optical signal delays exploiting the effect of Autler- 
Townes splitting in an inhomogeneously broadened quantum dot medium. The absorption and 
slow-down effects are compared for three schemes i.e. S, V and A, corresponding to different 
excitation configurations. Qualitative differences of the V-scheme compared to the S- and A-scheme 
are found, which show that features of Autler-Townes splitting are only revealed in the V-scheme. 
The underlying physical mechanisms causing this discrepancy are analyzed and discussed. Finally we 
compare field propagation calculations of the schemes showing significantly larger achievable signal 
delays for the V-scheme despite finite absorption of the coupling field. This opens the possibility 
for using waveguide structures for both coupling and probe fields, thus significantly increasing the 
achievable signal delays. 



I. INTRODUCTION 

A striking feature of a three level atom coherently 
driven by a strong coupling field and probed by a weak 
probe field, is a resonance splitting of the probed tran- 
sition, known as Autler-Townes splitting (ATS). The 
dressing of the states leads to a strong change of the 
probe absorption and phase; the latter corresponds to 
a change of the group velocity. A special case of ATS is 
electro magnetically induced transparency (EIT), 2 where 
the probe absorption vanishes within a narrow band due 
to destructive quantum interference among the different 
probe absorption paths. 

By exploiting EIT in a Na Bose-Einstein condensate, 
Hau et al. demonstrated in 1999 a tremendous reduction 
of the probe group velocity. Since then, astounding phys- 
ical phenomena based on EIT have been demonstrated, 
such as stored light, ' stationary light,' 1 single-photon 
slow-down and storage' and optical information process- 
ing with matter wave dynamics. A device enabling sig- 
nificant control of the velocity of light is of interest for ap- 
plications within optical communication and microwave 
photonics, see e.g. Ref. 9 and 10 and references therein. 
Compact and cheap devices, however, are required for 
practical applications, which makes semiconductor ma- 
terial a desirable material for slow light devices. 11 

Self-assembled quantum dots (QDs) are excellent 
candidates for realizing EIT in a solid material. 11 ' 1 1 QDs 
offer quantized energy levels with transition wavelengths 
that can be engineered. Furthermore, self-assembled 
QDs can be stacked to increase the QD density and em- 
bedded in waveguide structures to increase the laser field 
overlap. However, there are several issues that hinder 
the achievement of EIT in QDs: Firstly, dephasing rates 
in QDs are significantly higher than those found in atomic 
gases due to coupling to phonon modes. lo,1( ' While 
this problem may be dealt with by operating at cryo- 
genic temperatures, current fabrication techniques are 
inevitably associated with an inhomogeneous broadening 



(IHB) that, at cryogenic temperatures, is approximately 
three orders of magnitude larger than the homogeneous 
linewidth of the ground state transition of each individ- 
ual QD. IHB ~ J 1 is primarily associated with the size 
dispersion of the QDs 2 " and tends to destroy the observ- 
able EIT features when operating on an ensemble. 1 ' 5 ' 21 
As a result of these drawbacks, the associated slow light 
effects are typically mainly generated by the ATS rather 
than the more delicate interference effect of EIT. 

Lately, ATS in QDs have been reported using charged 
QDs, 22 as well as neutral QDs based on level schemes ex- 
ploiting either the fine structure splitting (FSS) 2 S ~ 25 or 
the biexciton transition 21 ' -2 ^ (see Fig. 8). One advantage 
of using two non-degenerate fine structure levels in a V- 
type configuration, is the ease of separating the pump 
and probe due to the linear orthogonal polarizations. 25 
Furthermore, as will be shown later, such a configuration 
could potentially be used in an inhomogeneously broad- 
ened QD ensemble. On the other hand, as pointed out 
in Ref. 29, a V-type configuration does not offer the ef- 
fect of EIT since the two absorption paths of the dressed 
states always interfere constructively. Rather, the re- 
sulting reduced absorption and group velocity is solely a 
result of ATS. Furthermore, the specific scheme based 
on fine structure split transitions, requires both fields 
to be polarized in the plane of growth '" (see Fig. 8). 
As a result, at least one of the two fields (coupling or 
probe) are required to propagate perpendicularly to the 
QD plane of growth. Letting the probe propagate per- 
pendicularly to the plane, offers the possibility of using a 
waveguide structure for the coupling field, thereby reduc- 
ing the required coupling power. However, this limits the 
interaction length of the probe field to only a few layers 
of QDs. Instead, if the waveguide structure is used for 
the probe, the coupling field must propagate perpendic- 
ular to the plane, and the beam profile needs to cover 
the full length of the waveguide. An unguided coupling 
field is impractical for integrated photonic devices, but 
equally important, the expanded mode area increases the 
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FIG. 1. The three EIT schemes 3 (left), V (middle) and A 
(right) as implemented for a quantum dot confinement poten- 
tial. Dashed blue and solid red arrows indicate the probe and 
coupling transition, respectively. 



power required for the coupling field. This problem is not 
present for the bi-exciton based scheme, where both fields 
are allowed to be either TE or TM polarized. However, 
as we shall see, this type of level configuration is partic- 
ularly sensitive towards IHB. Thus, to our knowledge, 
there are currently no experimental demonstrations of 
slow light based on ATS in a semiconductor waveguide 
structure. 

In this paper we propose a simple scheme that allows 
for ATS based slow light in an inhomogeneously broad- 
ened QD medium using a waveguide structure for both 
coupling and probe field. We start by introducing a sim- 
ple but general level structure for which the response of 
the ensemble of QDs is analyzed. The calculated slow- 
down and absorption effects of the three schemes are 
compared and discussed. The dependence on decay time 
constants for the V-scheme are further investigated in 
order to extract the underlying slow-down mechanisms. 
Finally, we present propagation simulations for the three 
schemes comparing the achievable signal delay for a given 
tolerated signal attenuation. 



II. QUANTUM DOT MODEL 

We consider QDs with discrete states, restricting our 
attention to two states in conduction band and a single 
state in the valence band. One transition is driven by 
a strong (coupling) laser field while a weak (probe) field 
probes a second transition. Three schemes S, V and A are 
considered that differ by the state used as a transit state, 
i.e. the state that is part of both the coupling and probe 
transition. The schemes are illustrated in Fig. 1. For 
ease of notation and comparison among schemes, we shall 
denote the three states in two ways: \a), \b) and |c) refer 
to the valence band state and the two conduction band 
states, respectively, while |1), |2) and |3) are denoted 
in accordance with the considered transitions. I.e. the 
probe always couples the states |1) and |3) while coupling 
field couples |2) and |3). 

Disregarding IHB the following features are noted: For 
a symmetric confinement potential, the V- and A-schcme 
involve a forbidden transition. I.e. in a symmetric po- 
tential, the transition |1) ■<-> |3) is forbidden due to the 
orthogonality of the envelope wavefunctions. U ' !J Thus, 



in order to exploit the V- and A-scheme, one needs to 
break the symmetry of the QD potential. This can be 
achieved, e.g. by applying an electric DC field across 
the QD. In reality, however, bandmixing and anisotropic 
strain effects relax these selection rules, thereby allow- 
ing for a finite coupling between the otherwise forbid- 
den transitions, 33,34 however, with a somewhat weaker 
dipole moment. For the A configuration this implies a 
reduction of the expected slow-down effect, while for the 
V-scheme a stronger coupling intensity is required for 
achieving the necessary ATS. As mentioned, it is only 
the H- and A-scheme that offer the possibility of EIT. 
As such, for a single QD, these schemes are expected to 
offer the most dramatic reduction of the group velocity 
and absorption. Especially the S-scheme seems as the 
preferred choice, since the probe transition, as opposed 
to the A-configuration, possesses a large dipole moment. 
Because the group index roughly scales with the dipole 
moment squared, a large dipole moment is desirable for 
slow light effects. 

The QD is modeled as an infinite disc shaped confine- 
ment potential with radius r and height In this case 
the undressed eigenenergies are found as: 13 




(1) 



where Wj is the effective electron or hole mass and Z n \ 
denotes the n'th rooth of the Vth Bessel function, i.e. 
Ji(Zni) = 0. Since the QD sizes are dispersed in the 
ensemble, we explicitly write the radius, r, and height, 
<;, dependence of the eigenenergy in Eq. (1). Assuming 
all QDs to have a constant aspect ration r/ — s/r, an 
infinitesimal volume change of the QD with a change of 
radius by Sr results in a shift of the eigenenergy of 
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Due to the heavier effective mass of the hole states, we 
consider only energy shifts of the electron states. Setting 
i] = 3/10 18 the shift of the first two electron states are 
found as 



feioi 



fen 



-115.4 £r 



-62.17 Sr. 



(3) 
(4) 



Linearizing the problem is justified since the inter-band 
transition energy is typically ~ 1 eV, and the inhomoge- 
neous broadening is on the order of a few tens of meV. 18 
The shift of eigenenergies leads to a shift of the opti- 
cal transitions that we shall refer to as a spectral shift 
and is illustrated in Fig. 2. For a given subensemble 
of QDs with fixed energies, we denote the spectral shift 
of the probe (coupling) transition p (Aj/ lc ) and de- 
fine it relative to the mean probe (coupling) transition 
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FIG. 2. Density of resonant transitions as a function of fre- 
quency for the three considered transitions indicated by the 
level diagram insets. The IHB FWHM is in this work referred 
to the ground state transition (middle peak). 



energy u)i 3 of the QD ensemble, i.e. Aj/j )P = W13 — Q13 
(A,/, lC = w 23 — W23). Relating the spectral shift of the 
probe transition, A,-^ iP , to the spectral shift of the cou- 



pling transitions A,-^ iC we find for the three schemes: 

S: A ih , c = «&*^ml A ih)P = K A ih , p « 0.077 A ife , p (5a) 
V: A lh . c = ff^A iftiP = nA lh . p a 1.08A i/l)P (5b) 
A: A i/t , c = 3(£l ^ 10l) A^, p = kA iKp a 0.072A ih , p (5c) 

where « = A ih , c /A lhtP . 

The absorption coefficient and group index are re- 
lated to the imaginary and real part of the electric sus- 
ceptibility Xp associated with the probe transition. \p 
is calculated from the density matrix equations assum- 
ing monochromatic classical electric probe and coupling 
fields and phenomenologically adding population and po- 
larization decay rates." For one sub-ensemble of QDs with 
fixed energy levels (size), \p calculated to first order in 
the probe field evaluates to: 



V: Xp (A p ,A c ) 
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where V is the QD volume, T is the confinement fac- 
tor and £0 is the vacuum permittivity. A p and A c 
are the effective probe and coupling laser detuning, i.e. 
A„ = A.„ 

and A c = uj c — |w 3 — cj 2 |. A two-photon 
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A p ± A c , fiis is the dipole 



detuning is defined as S± 
moment of the probe transition and the Rabi frequency, 
il c of the coupling field is defined as il c = H23 ■ ^c/(2S) 
where £ c is the coupling electric field amplitude. Popu- 
lation and dephasing (polarization) decay rates between 
the states i and j are denoted Tij and 7^ , respectively. In 
the expression for the V scheme we furthermore defined 
Ti = r i3 - r i2 and T 2 = T 23 + r i2 , i.e. the total pop- 
ulation loss rate of state 1 and 2, respectively. Finally, 
we defined £ = 2Fi 3 r 2 (A 2 . + 723)- Comparing equations 
(6b), (6c) and (6a), it is seen that the susceptibility only 
depends on the homogeneous linewidths for the S- and 
A-scheme, whereas for the V-scheme the population de- 
cay rates also become important. As we shall see later, 
this puts restrictions on the population decay rates in 
order to achieve ATS using the V-scheme. Denoting the 
inhomogeneous distribution function of the probe transi- 
tion by /, the mean electric susceptibility is obtained by 



averaging over the QD ensemble: 

00 

(x(A p ))= J f(A ih ) x {A p -A ih ,A c -KA ih )dA ih , (7) 



where the subindex p was omitted in A^. 

To the author's knowledge, no publications are avail- 
able, giving the measured population decay and dephas- 
ing rates on the same QD sample for transitions related 
to those depicted in Fig. 1. Thus, the population de- 
cay rates are set similar to those reported by Heitz et 
al.. 18 We note though, that in this paper they report 
on QDs with an exceptionally low intra-band popula- 
tion decay rate. It is demonstrated, that the geometrical 
shape of the QDs can be tailored to either enhance 35 or 
reduce the exciton-LO-phonon coupling. In the latter 
case, truncated pyramidal shaped QDs where used, giv- 
ing measured intra-band relaxation rates ~15 times lower 
than the radiative population decay rate of the ground 
state transition. Later, we shall discuss the consequences 
of using typical QDs with a much faster intra band de- 
cay rate. For the |c) —> \a) decay, we tentatively set 
r ac = r a b. We note, though, that for the considered 
case of a closed 3-level system, the transition rate r ac 
is governed by the dipole moment squared, cf. Fermi's 
golden rule. M ' Since fi ac is set very low, a lower value of 
r ac relative to T a b is predicted by Fermi's golden rule. 
In reality however, carrier-carrier scattering and carrier- 
phonon interaction 1 ' lead to larger decay rates, i.e. the 
system is not truly a closed 3-level system. 

Concerning the dephasing rates, earlier measurements 
show near lifetime limited decay rates of the ground state 
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Parameter |a) -o |&) \b) o |c) |a) -o |c) 

»7 3/10 

r 6-10 -3 

V 1200 nm 3 

(Tift 10 meV 

hutij 0.999 eV 61 meV 1.06 eV 

fHj 0.7 e-nm 4.7 e-nm 0.10 e-nm 

Tij 2.6 A*eV 0.16 ^eV 2.6 jieV 



TABLE I. Quantum dot parameters used in the calculations 
unless otherwise stated. r\ is the average height-radius ratio, V 
is the confinement factor, V is the average QD volume and 
is the FWHM of the spectral distribution function /. hujij, 
Hij and Tij are the energy difference, the dipole moment and 
the population decay rate of the transition between \i) and 
The homogenous linewidths 7^ are all assumed lifetime 
limited. 

transition near K temperature. 15 ' 38 ' 39 Thus, for the cal- 
culations we tentatively set the dephasing rates of all 
transitions to be lifetime limited unless otherwise stated. 
The parameter values are given in table I. 

III. SLOWDOWN AND ABSORPTION 
CALCULATIONS 

A. Scheme comparison 




Pump Intensity [W/cm 2 ] 



FIG. 3. (color online) Top: Normalized absorption as a func- 
tion of coupling intensity for the three schemes H (red), V 
(black) and A (blue), with (solid) and without (dashed) IHB. 
The absorption is normalized by the absorption without an 
applied coupling field. Note that the curves for S and A over- 
lap. Bottom: corresponding slowdown factor S. Material 
parameters are presented in table I. 



As a figure of merit describing the degree of group 
velocity control, we define the slow-down factor as S = 
Ug/nbg, where n^g is the refractive index of the surround- 
ing lossless background material while n g is the calcu- 
lated group velocity including the QD material and given 
as 

n g = Rc (yi + Xbg + Xp) + ^—q^ . 

(8) 

where Xbg is the electric susceptibility of the background 
material. In Fig. 3 the probe absorption (top) nor- 
malized by the absorption without a coupling field and 
the probe slow-down factor (bottom) is seen plotted as 
a function of coupling intensity for the three different 
schemes using (6b)- (6a). Both probe and coupling de- 
tuning are set on resonance with the average transition 
frequencies W13 and H123, respectively. Calculations with 
(solid) and without (dashed) IHB are presented. Com- 
paring the three schemes without IHB, the maximum 
slow-down of the V-scheme requires a stronger coupling 
intensity compared to the other two scheme. This is 
a result of the weak dipole moment between \a) and 
\c). Thus, for the V-scheme, a fairly high intensity is 
required in order to realize the required ATS. The A- 
scheme, shows a somewhat lower slow-down factor as a 
result of the low dipole moment of its probe transition. 
From the absorption plot, we note that the S- and A- 
schemes overlap due to the identical dipole moment of 



the coupling transitions and the similar dephasing rates 
712 and 7 13 , see equation (6b) and (6c). When compar- 
ing the calculated slow-down factors including IHB, it 
is seen that the V-scheme is far superior. For the two 
other schemes, the inhomogeneous broadening is seen to 
completely destroy any slow-down effect, whereas the V 
scheme maintains a slow-down factor of more than 10 2 . 
Also, the required coupling power for optimal slow-down 
is smaller by roughly 3 orders of magnitude, relative to 
the other two schemes. Finally, the V-scheme shows an 
almost similar absorption as a function of coupling inten- 
sity when including IHB. This certainly is not the case 
for the S- and A-scheme, where the coupling field needs 
to be increased by more than 8 orders of magnitude to 
achieve the same reduction of absorption. 

To understand these notable differences we need to dis- 
cuss the details of the ATS mechanism. The presence of 
a strong coupling laser turns the bare eigenstates \2) and 
1 3) into the dressed states |+) and |— ) that appear as two 
absorption resonances for the probe. In a rotating frame 
the interaction Hamiltonian can be written: - 
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where A x = A p - A ih (A 1 = -A p + A ih ) for the V (3 
and A) scheme and A2 — A c — KAih (A2 = -Ac + kA^) 
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for the V and 3 (A) scheme. Letting il p — > we find the 
eigenvalues as 

A 1 = -A 1 , A ± = ^-A 2 ± v /4fi2 + A^ (10) 

where the X± are the shifts of the dressed resonances of 
|3). Thus, the resonance condition for the probe transi- 
tion is shifted accordingly. For a class of QDs spectrally 
shifted by A^, this becomes 

5, A: A p = A ih + A± V: A p = A ih - A± (11) 

Suppose the probe laser is on resonance with QDs at the 
center of the distribution function i.e. A p = A^ = 0. 
Since the coupling field has shifted the resonance by A±, 
QDs at the center of the distribution function do not 
contribute with absorption, assuming the coupling Rabi 
frequency larger than the linewidth of the probe transi- 
tions. However, for A^ ^ this may not be the case. 
Rewriting (11) for A p = 0we find: 

([2± K ] 2 - K 2 )A^4^, (12) 

where the sum is for the 5-scheme and the difference is 
for the V- and A-scheme. A real solution to (12) for A^ 
exists when 

H:k>-1 and V,A: k<1 (13) 

Comparing with the calculated values of k in (5), using 
the discpotential, we see that (13) is only fulfilled for 
the 3- and A-scheme. Physically, this means that for 
the V-scheme, none of the QDs in the inhomogeneously 
broadened ensemble have one of the two dressed states 
on resonance with the probe, as opposed to the other 
two schemes. Importantly, this implies that the probe in 
the V-scheme does not experience significant absorption 
from any of the QDs in the ensemble. To this we may 
add that we expect this observation holds true even when 
using a more realistic QD model, e.g. by using a more re- 
alistic QD potential as well as including the energy shifts 
of the holes. Generally, one finds that k > 1, k > 1 and 
< k < 1 for the V- and A-scheme respectively. This 
difference is clearly seen in Fig. 4 where the real and 
imaginary parts of the electric susceptibility Xp a re plot- 
ted as a function of A p and A^. Left, middle, and right 
column are the electric susceptibility of the V- and 
A-schemes, respectively, with the top and bottom plots 
showing the real and imaginary part of the susceptibility 
X P normalized by [eo^(r/i 2 ) -1 ] . Without loss of gener- 
ality, we have in this plot, for the purpose of illustration, 
used a height-radius aspect ratio r\ — 1, although this 
ratio is typically found smaller. 18 The plot clearly shows 
two resonances crossing the region where A p = for the 
A and S scheme. Thus, the average susceptibility, plot- 
ted in the top row of Fig. 4 using equation (7), reveals 
that the transparency region is destroyed and slow light 
effects vanish. 




V r i3 V F 13 

FIG. 4. Electric susceptibility (in units of r fif 3 / '(V Eoh)) for 
QDs with rj = 1. Left middle and right column is for the S, 
V and A, respectively. Top: Average real (dashed) and imag- 
inary (solid) part of the susceptibility as a function of the 
normalized probe detuning. Middle: Real part of the electric 
susceptibility as a function of the normalized probe detun- 
ing, A p , and spectral shift A^. Bottom: Imaginary part 
of the electric susceptibility as a function of the normalized 
probe detuning, A p , and spectral shift A^. The asymptote of 
the primary and secondary resonance is indicated by a dash- 
dotted and dashed line, respectively. 



An illustrative way of exposing the difference between 
the schemes is by considering the resonance condition 
(11) for QDs largely detuned from the center of fih, i.e. 
\Aih\ 3> fi c . Using the limit of (10) for |A 2 | — > oo, we 
find the two approximate resonance conditions: 

I: A th = A p and II: A ih = (1 ± n)' 1 A p , (14) 

where the sum is for the S-scheme while the difference is 
for V- and A-schemes. We denote the two resonances the 
primary and secondary resonance. ~ ' The two asymptotes 
are illustrated as a solid and dashed line in Fig. 4. By 
virtue of the values of k in (5), the asymptotes of A^ 
vs. Ap have a positive slope for the 2- and A-scheme. 
From Fig. 4, this clearly implies that for these schemes, 
the resonance crosses the region where A p = 0. The 
V-scheme, however, exhibits a negative slope, and as a 
result, no dressed state is resonant with the probe field. 

Another notable difference between the V-scheme com- 
pared to the two other is that only near A^ « do two 
resonances occur, i.e. the secondary resonance vanishes. 
In the bare state picture, the primary resonance is inter- 
preted as the absorption from |1) — > |3) while the sec- 
ondary resonance is, for S and A, the sum of absorption 
pathways |1) (|3) ->• |2) ->• |3)) n where n is the number 
of Rabi floppings on the coupling transition. For the V- 
scheme the pathways are (|3) — > |2) — > |3)) n — > |1), where 
the coupling field is seen to excite carriers. Thus, for 
those QDs where the coupling laser is far off resonance, 
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FIG. 5. Top: absorption as a function of coupling field inten- 
sity for the three cases: (a, solid) lifetime limited dephasing 
rates with no intraband population decay rate, (b, dashed) 
lifetime limited dephasing rates with a large intraband pop- 
ulation decay rate and (c, dotted) being equivalent to (a) 
except the intraband dephasing is set large. Bottom: corre- 
sponding slow-down factors. 



i.e. | Ay, | 3> 0, carriers from |3) are no longer excited to 
1 2). As a result, the absorption pathways associated with 
the secondary resonance vanishes in the V-scheme. 

B. V-scheme parameter dependence 

As mentioned, we have chosen a very low intraband 
decay rate ~ 8 ns. Although such QDs have been 
fabricated, 1 "' quantum dots typically show intraband 
relaxation on a time scale of ~ 10 ps. In Fig. 5 we 
compare the absorption and slow-down of the V-scheme 
in three cases corresponding to different decay time con- 
stants: 

(a) Lifetime limited polarization decay with low intra- 
band population decay rate; ri2/ri3 = 0. 

(b) Rates equivalent to a) except the intraband polar- 
ization decay 712 is set a thousand times larger 
than the lifetime limited value, i.e. 712 = 1000 • 

5(ri3 + r 23 + r 12 ). 

(c) Lifetime limited polarization decay and a large in- 
traband population decay rate; T12/T13 = 25, cor- 
responding to a relaxation time of 10 ps. 

Comparing the two cases (a) and (b) in Fig. 5 , we 
notice an almost similar slow-down factor for large cou- 
pling intensities, but for small intensities the slow-down 
for case (a) is approximately 2.5 larger than that of (b). 
This gives insight about the origin of the slow-down ef- 
fect in the V-scheme: Since the polarization rate in (b) is 




FIG. 6. Plot of real (top) and imaginary (bottom) part of 
the electric susceptibility \ p (in units of r/if 3 / '(Veoh)) as a 
function of the scaled probe detuning A p /Fi3 and intraband 
decay rate Fi2/ri3. 



set very large, the observed slow-down is not due to AT 
splitting since this requires a coupling Rabi frequency 
effectively similar to 712. Rather, the slow-down origi- 
nates from spectral hole burning (SHB) caused by the 
coupling field. The hole burning causes the probe to ex- 
perience a reduced absorption near its resonance which 
leads to a positive slope of the refractive index. This is 
confirmed by the variation of the absorption, which is 
seen to decrease to half its value for increasing coupling 
intensity. When increasing the Rabi frequency beyond 
712 (~ 10 7 W/cm 2 ), the absorption decreases further to- 
wards zero due to ATS. Thus, the observed slow-down 
effect of the V-scheme for lifetime limited decay, (a), is 
partly due to AT splitting and partly due to spectral 
hole burning. Notably, this is observed in Fig. 3, at low 
coupling intensities. Here, the slow-down, for decreas- 
ing coupling intensity, is seen to decrease more gradually 
than for the other schemes, due to the SHB origin of the 
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FIG. 7. (color online) Calculated electric susceptibility using 
rj = 1. Top: Average real (dashed) and imaginary (solid) part 
of the electric susceptibility (normalized by Tfii 3 /(Veoh)) for 
the alternative V scheme (red) and V-scheme from Fig. 1 
(grey). Bottom: Associated imaginary part of the electric 
susceptibility for the alternative V-scheme as a function of 
scaled probe detuning A p /Ti3 and spectral shift Aih/Ai3. 
The inset illustrates the excitation configuration of the alter- 
native scheme. . 



slow-down. 

Considering the case of a large intraband population 
decay rate, (c), we notice that a larger coupling inten- 
sity is required at the maximum slow-down. This is a 
result of a larger value of 712 given that is a sum of all 
three population decay rates. As a result, a larger Rabi 
frequency is required to resolve the AT splitting. For a 
large intraband decay rate, the system acts as an ampli- 
fier with electrons pumped by the coupling field from |3) 
to 1 2) and quickly decaying to |1). Thus, for large cou- 
pling intensities (~ 10 4 W/cm 2 ), the probe experiences 
gain (negative absorption) resulting in a slow-down ef- 
fect. For the largest intensities, > 10 4 W/cm 2 , the ab- 
sorption is seen to return to transparency for increasing 
coupling intensity as a result of the ATS. This is il- 
lustrated in Fig. 6 where the real (left) and imaginary 
(right) part of the susceptibility is plotted as a function 
of the probe detuning, for various intraband scattering 
rates r 12 . For T 12 = 0, the absorption spectrum shows 
ATS as a result of the strong coupling field. However, as 
Ti2 is gradually increased, the absorption becomes nega- 
tive (i.e. amplification). Considering the real part of \ v 
(left), it is seen that the slope remains positive, giving a 
reduction of the group velocity even for the highest val- 
ues of r 2 3. It is interesting to note though, that the two 
effects support each other. Thus, for intermediate val- 
ues of r 12 an effective slow-down may be achieved, while 
minimizing the absorption due to the small fraction of 



carriers being pumped to the probe transition. 



C. Alternative V schemes 

The V-scheme may also be considered in an alternative 
configuration where the probe and coupling fields have 
interchanged transitions as illustrated in the inset of 
Fig. 7. We shall denote this the alternative V-scheme. 
However, there are several reasons for choosing the 
standard V scheme illustrated in Fig. 1 rather than the 
alternative configuration. Firstly, the expected weak 
dipole moment of the "forbidden" transition would 
result in a small slow-down of the probe, since it scales 
with the dipole moment square. Also, a finite intraband 
population decay rate would lead to an increased probe 
absorption as well as a weakened slow-down. More 
importantly, however, is the associated value of k. For 
the alternative V-scheme we have k 0.43. Thus, 
real solutions to (11) exist, i.e. the slope (1 — k) _1 
is positive (see equation (14)) leading to an increased 
absorption of the probe field. This is illustrated in 
Fig. 7 (bottom) where the probe absorption is plotted 
as a function of probe detuning A p and spectral shift 
Aih- The two dressed absorption resonances are seen 
to overlap the region with A p = 0, and the average 
spectral dip of the probe absorption near resonance 
thus becomes less pronounced, resulting in a lower 
slow-down factor. This is clearly seen in Fig. 7 (top) 
where the absorption integrated over is plotted 

as a function of A p for the two types of V schemes. 
This clearly demonstrates, that it is not only the 
SHB slow-down effect that makes the V-scheme superior 
to the other schemes, but also the value of n is important. 



D. Fine structure splitting schemes 

As mentioned, ATS in QDs has previously been 
demonstrated by exploiting the FSS in QDs. The FSS 
arises from the electron-hole exchange interaction that 
couples the electron spin and hole angular momentum 
creating two dark excitonic states with total angular mo- 
mentum projection j z = z + j e z = ±2 and two bright 
states with j z = jh.z+je.z = ±1. 12 ' 30 ' 41 The electron- hole 
coupling leads to a splitting between the dark and bright 
states, and we shall therefore neglect the former. In a QD 
with rotational symmetry in the plane of QD growth, the 
two bright states are degenerate. However, typically QDs 
have elongated potentials along the crystal axis [110] and 
[110], either due to the structural formation of material 
or due to the strain induced piezoelectric fields. 23 ' 30 ' 41 
When the rotational symmetry of the QD confinement 
potential is broken, the angular momentum is no longer a 
good quantum number and the bright eigenstates evolve 
into a symmetric and antisymmetric linear combinations 
of the j z = ±1 basis states, with opposite energy shifts 



[110] 





FIG. 8. a) Asymmetric QD potential with symmetry axis 
(dashed lines) along the crystal planes [110] and [110]. b) Fine 
structure splitting based schemes. The two exciton states 
\X\\) and \X_i_) are coupled with linearly polarized light par- 
allel to the symmetry axis. 



I^fss- 12 We shall denote these eigenstates \X±) and 
since they are coupled by light that is linearly polarized 
perpendicular and parallel to the two axis of reflection 
symmetry of the QD, see Fig. 8a. Considering the biexci- 
ton state \XX), where two excitons are excited, we note 
that this state is not split by the exchange interaction 
since the net angular momentum of both excitons is 0.' i(l 
Together with the crystal ground state |0), these states 
may be optically driven at the transitions that make up 
either the S- V- or A-scheme, by suitable choice of wave- 
length and polarization of the coupling and probe field 
(see Fig. 8). We shall therefore extend our model to 
these types of schemes. 

Similar to equations (5), relating the spectral shift of 
the coupling transition to the probe transition we find 
for the V-scheme: 



V: 



S (e x ± ^£fss) 
8 (ex T ^fss) 



±ih 7 p 



l ± 



FSS 



8 (e x T §£fss) 



(15) 
(16) 
(17) 



where ex denotes the exciton energy without the FSS 
and in the last equation it was assumed that ex S> eFSS ■ 
The FSS, eFSS) has previously been reported to vary with 
the size of the quantum dots, 5(1 i.e. the FSS decreases 
with increasing groundstate transition energy. However, 
the magnitude compared to the IHB is rather low, e.g. 
epss w &s found to change approximate 0.5 meV for a 
shift of the ground state transition of roughly 300 meV. 
Thus, the last term in (17) is neglected leading to a value 
of K « 1. 

Analogously we find that k~1 for the A-scheme. 
For the 2-scheme we get: 



5e 



xx 



Se x 



= 1 



5 (e X x - ex 
Sex 



A 



ih,p 



(18) 
(19) 



Generally, the biexciton energy, exx, is different from 
that of the exciton, ex, due to mixture of Coulomb inter- 
action and correlation effects that depends on the quan- 
tum dot size, shape and strain. Thus, a change of QD 
size does not in general lead to identical shifts of the exci- 
ton and biexciton transition energies, i.e. Sexx/Sex 7^ 0. 
In a series of measurements by S. Rodt et al, they mea- 
sured the energy difference exx — ex as a function of 
ex- Applying a linear fit to these measurement we find 
S (exx - ex) /Se x « 0.05. 

Comparing the values found for k with the conditions 
in (13), it is immediately seen that the S-scheme remains 
unsuitable for achieving slow light since some QDs con- 
tribute to the probe absorption. However, for both the 
A- and V-schemes a spectral region near A p = remains 
transparent since no dressed state is on resonance with 
the probe. As mentioned previously, a general problem 
of exploiting a FSS scheme is that it requires the cou- 
pling field to propagate perpendicular to the plane of 
growth. With the A-scheme we are furthermore left with 
the problem of population decay from the exciton state, 
\X), to the crystal ground state |0). Thus, a prepulse 
is required to optically pump a population into \X). As 
pointed out in Ref. 43, the prepump needs to arrive suf- 
ficiently long time prior to applying the probe field to 
avoid dipole interference induced by the prepump. This 
time is governed by the dephasing rate of the |0) — \X) 
transition. On the other hand, the prepump should not 
arrive earlier than the relaxation time between \X±) and 
\X\\), to avoid population in both of the exciton states. 
Thus, a proper calculation of the electric susceptibility 
requires solving the density matrix and wave equation 
in the time domain. However, restricting our attention 
to the impact of inhomogeneous broadening, we shall for 
the A-scheme neglect population decay from \X) to |0). 

The calculated imaginary part of the electric suscep- 
tibility is plotted in Fig. 9 for the V-scheme (left) and 
A-scheme (right) as a function of the probe detuning and 
spectral shift. In both cases we find k — 1, but the sec- 
ondary resonance remains absent for the V-scheme due 
to carrier excitations by the pump, as discussed in sec- 
tion III A. When integrating over all spectral shifts, this 
leads to significantly different average susceptibilities for 
the two schemes (Fig. 9 top). Where the V-scheme looks 
similar to the calculation based on the conduction band 
states in section III A, the A-scheme shows a very sharp 
and narrow dip in the absorption and consequently a 
steep slope of the real part of the susceptibility. As a 
result, the FSS based A-scheme could potentially be used 
to realize very low group velocities. 



IV. PROPAGATION EFFECTS 

Since the V-scheme turns out to be robust against IHB, 
the use of a waveguide structure where the electric fields 
interacts with an ensemble of QDs, is possible. The pur- 
pose is two fold: the confinement of the electric fields 
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FIG. 9. (Color online) Electric susceptibility (in units of 
T [i\ 3 / (V Eoh)) of the V (left) and A (right) scheme based on 
the FSS configuration. Bottom: Imaginary part of the electric 
susceptibility as a function of the normalized probe detuning, 
A p , and spectral shift Aj^. The asymptote of the primary and 
secondary resonance is indicated by a dash-dotted and dashed 
line, respectively. Top: Corresponding calculated average real 
(dashed) and imaginary part (solid) part of the electric sus- 
ceptibility. Key parameters are T p = r c and Y uc = 0, where 
T p , r c and T uc denote the population decay rate of the probe, 
coupling and uncoupled transition, respectively. Dephasing 
rates were assumed lifetime limited. 



indicate the z-dependence of Xp due to the absorption 
of the coupling field in the V scheme. Equation (20) is 
formally solved to give: 




{x P {z,u))dz 



(21) 



Noting that (21) takes the form £{z) — £(0) exp(ifc), 
where k is the mean propagation wavevektor over the 
propagation length z, we can extract the average group 
index (n g ) and absorption coefficient (a) using the rela- 
tions rig/c — dRe(k)/duj\ LUp and a = Im(fc). Denoting 
the real and imaginary part of \ as \' an d x" ■> respec- 
tively, we find the average group index as: 

KM) = — + 



' I (x qd ( z ,u)) +uj H — j dz, (22) 



1 1 

z 2no 

while the average absorption coefficient is found as: 



1 u p 

Z 271kg C 



(Xp(z,uj))dz. 



(23) 



serves to decrease the nessesary pump power, and fur- 
thermore the waveguide provides an efficient way of in- 
creasing the interaction length, thereby increasing the 
potential signal delay. However, as opposed to the two 
other schemes, the coupling field is subject to carrier ab- 
sorption. An upper limit on the achievable delay is there- 
fore expected. I.e., for a long waveguide the coupling field 
is initially very strong and generates a large ATS in the 
first part of the waveguide. The probe absorption in this 
section is very low and so is the slow-down effect. At the 
end of the device, on the other hand, the coupling field is 
absorbed such that the limited ATS causes a high degree 
of absorption, while the slow-down effect on the probe 
has vanished. Since this problem is avoided with the El- 
and A-scheme, a comparison of the delay and transmis- 
sion of the three schemes is reasonable. 

The propagation effects of both the coupling and probe 
fields were analyzed using a one-dimensional propagation 
equation, assuming a single transversal mode. 

The propagation of the probe field was calculated 
from the wave equation in the slowly varying enve- 
lope approximation. 5 In the frequency domain the wave 
equation for the probe field is given as: 



d£(z,n) 
dz 



n 



2n hn , 



( Xp (z)))£(z,n), (20) 



where z is the propagation distance, £ is the slowly vary- 
ing electric field amplitude, u g is the group velocity of 



Notice that equation (22) and (23) describe the average 
group index and absorption, respectively, over the prop- 
agation length z. Using (22) the aquired probe delay is 
calculated as 



At = - {{n g [z)) - n bg ) 



(24) 



For ease of calculations, the propagation of the coupling 
field, in the V-scheme, was decoupled from that of the 
probe field by approximating the pump transition to be 
a pure 2-level problem, with two decay channels: One di- 
rectly from 1 2) to |3) and another via the state |1), where 
the population decay rate for the latter channel is given 

by (r 
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the background material and f2 



Dp. We explicitly 



This approximation is valid for a 
weak probe and Ti2 <C T23, meaning that the population 
density of 1 1) is vanishing. 

In Fig. 10 the delay and transmission are plotted as 
a function of propagation length z and injected coupling 
intensity for all three schemes 5, V and A. The surface 
plot shows the delay in units of ns, while the lines on 
top show iso-curves of the calculated transmission in dB. 
Noting the different scalings of the associated colors bar, 
it is seen that for the 5 and A scheme, hardly any posi- 
tive delay is possible as a result of the large IHB. For the 
V scheme, however, a considerable delay, on the order of 
nanoseconds, is obtained. Furthermore, it is seen that the 
required coupling intensity for maximum delay, is consid- 
erably larger for both the S and A scheme compared to 
the V scheme. In fact, such high coupling intensities in- 
evitably requires a strong field confinement such as in 
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Coupling intensity [W/cm ] 

FIG. 10. (Color online) Calculated delay (surface plot) and 
transmission (isocurves) as a function of injected coupling in- 
tensity and propagation length z for the A- (top), V- (middle) 
and H-scheme (bottom). The transmission curves are labelled 
with the corresponding transmission in dB. Note from the 
colorbar the different scales of the delay. Decay and material 
parameters are chosen as in table I. 




Coupling intensity [W/cm 



FIG. 11. Delay as a function of the injected coupling inten- 
sity at a fixed transmission of —10 dB for the V- (solid) and 
H-scheme (dashed). Corresponding propagation lengths are 
indicated by arrows for the V- scheme. Inset shows a zoom-in 
of the calculated delay for the H-scheme. Decay and material 
parameters are chosen as in table I. 



a waveguide. However, for these schemes, the coupling 
transition is typically around 10 /im whereas the probe 
transition is roughly 1 /im. Fabricating a waveguide for 
both wavelengths is a challenging task. 

Considering the calculated transmission coefficients, 
the A-scheme shows a considerably larger transmission 
compared to the other two schemes as a result of the 
small dipole moment of the probe transition. However, 
this also affects the achievable delay, being very small. 

Inevitably what sets the limit of the achievable delay, 
is the tolerated absorption of the probe signal. A plot of 
the delay as a function of coupling intensity along the -10 
dB transmission contour line is seen in Fig. 11 compar- 
ing the V and S schemes. From the inset it is seen that 
the delay of the S scheme continues to grow for increas- 
ing coupling intensity. This is due to the fact that the 
coupling field does not experience any absorption. Thus, 
increasing the coupling intensity results in a decreased 
delay that is compensated by increasing the propagation 
length. However, group velocity dispersion and wave- 
guide losses, that are not included in the model, set an 
upper limit on the device length. 1 ' ' r As opposed to the 
S- and A-scheme the delay obtained using the V-scheme 
is seen to approach an upper limit. Thus, the achieved 
delay is almost constant for intensities above 500 W/cm 2 
only changing by less than 6 %. Considering the device 
length, indicated by arrows, this implies that very lit- 
tle is gained by making the device longer than ~ 1 mm. 
Furthermore, it was found that this upper limit roughly 
scales linearly with the tolerated transmission (in dB). 
For a Rabi frequency larger than the dephasing rate we 
have the approximate relations At oc zn g (z) oc zf2 _1 . 



The AT splitting and therefore also the absorption a 
is governed by the Rabi frequency of the coupling field. 
For increasing Rabi frequency, the absorption converges 
towards zero. Since the maximum delay was found to 
scale with the logarithm of the tolerated transmission, 
this implies that (a) oc £2 . Thus, for a given toler- 
ated transmission, the field strength of the coupling field 
should have a magnitude where the average absorption 
coefficient is roughly inversely proportional to the Rabi 
frequency in order to achieve the largest delay. 

The above calculations are drastically simplified by 
only considering propagation of monochromatic waves. 
As such, the calculated delay serves as an upper limit, 
since group velocity dispersion as well as absorption dis- 
persion may restrict the device length. 4 '' 48 



V. CONCLUSION 

Slowdown effects in an inhomogeneously broadened 
quantum dot medium were analyzed for an ensemble of 
three level QDs strongly driven by a coupling field and 
probed by a weak field. 

Using a simple but general model for the quantum 
dot energy levels, the slow-down and absorption of 
three generic schemes, 5, V and A, were compared. 
When neglecting inhomogeneous broadening, the IB- 
scheme seemed advantageous both in terms of slow-down 
and absorption compared to the other schemes. However, 
when including inhomogeneous broadening all slow-down 
effects vanished for the EE- and A-schcmc, while the V- 
scheme maintained a slow-down of more than 10 2 . 

Comparing the dressed state eigenenergies it was 
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shown that no quantum dot in the inhomogeneously 
broadened ensemble would have a dressed state on reso- 
nance with the probe in the V-scheme. Thus, as opposed 
to the other schemes, a transparency region remains near 
zero detuning of the probe that results in reduction of 
the group velocity. By analyzing the V-schemc for differ- 
ent intraband relaxation rates, it was further shown that 
origin of the slow-down effect was an admixture of the 
Autler-Townes splitting and spectral hole burning caused 
by the coupling field. 

Finally, based on field propagation calculations the 
transmission and signal delay of the three schemes were 
compared. Despite the finite absorption of the coupling 
field, the V-scheme showed significantly larger achievable 
delays. For a fixed transmission of -10 dB the delay was 
seen to converge towards a constant value (~ 1 ns) for an 



increasing coupling intensity. As a result, the delay was 
almost constant for propagation lengths beyond 1 mm. 

The theory and simulations presented in this paper 
indicate the possibility, by proper tailoring of the level 
scheme and level lifetime, of experimentally achieving op- 
tically controlled slow light in an inhomogeneously broad- 
ened quantum dot medium using a waveguide structure 
for confinement of both probe and coupling field. 
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